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Abstract A study of W Z production in proton-proton 
collisions at y/s = 7 TeV is presented using data cor- 
responding to an integrated luminosity of 4.6 fb -1 col- 
lected with the ATLAS detector at the Large Hadron 
Collider in 2011. In total, 317 candidates, with a back- 
ground expectation of 68 ± 10 events, are observed in 
double-lcptonic decay final states with electrons, muons 
and missing transverse momentum. The total cross- 
section is determined to be a\^ z — 19.0l 1 ; 3 (stat.) ± 
0.9(syst.) ± 0.4(lumi.) pb, consistent with the Standard 
Model expectation of 17.6^J'Jpb. Limits on anoma- 
lous triple gauge boson couplings are derived using the 
transverse momentum spectrum of Z bosons in the se- 
lected events. The cross-section is also presented as a 
function of Z boson transverse momentum and diboson 
invariant mass. 



1 Introduction 

The underlying structure of electroweak interactions in 
the Standard Model (SM) is the non-abelian SU(2) L x 
U(l)y gauge group. This model has been very success- 
ful in describing measurements to date. Properties of 
electroweak gauge bosons such as their masses and cou- 
plings to fermions have been precisely measured at LEP, 
the Tevatron and SLD [T]. However, triple gauge boson 
couplings (TGCs) predicted by this theory have not yet 
been determined with a similar precision. 

In the SM, the TGC vertex is completely deter- 
mined by the electroweak gauge structure and so a pre- 
cise measurement of this vertex, for example through 
the analysis of diboson production at the Large Hadron 
Collider (LHC), tests the gauge symmetry and probes 
for possible new phenomena involving gauge bosons. 
Anomalous TGCs, deviating from gauge constraints, 
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Fig. 1 The SM tree-level Feynman diagrams for W^Z pro- 
duction through the s-, t-, and ti-channel exchanges in qq' 
interactions at hadron colliders. 

may enhance the W ± Z production cross-section at high 
diboson invariant masses. The cross-section can also be 
enhanced by the production of new particles decaying 
into W ± Z pairs, such as those predicted in supersym- 
metric models with an extended Higgs sector and mod- 
els with extra vector bosons [2]. 

At the LHC, W ± Z diboson production arises pre- 
dominantly from quark-antiquark initial states at lead- 
ing order (LO) and quark-gluon initial states at next-to- 
leading order (NLO) (3j. Figure [1] shows the LO Feyn- 
man diagrams for W^Z production from qq' initial 
states. Only the s-channel diagram has a TGC vertex 
and is hence the only channel to contribute to potential 
anomalous coupling behaviour of gauge bosons. 

In proton-proton (pp) collisions at a centre-of-mass 
energy yfs — 7 TeV, the SM cross-section for W^Z 
production is predicted at NLO to be 17.61^'J pb. This 
has been computed for 66 < m« < 116 GeV, where mu 
is the invariant mass of the dilepton system from the Z 
boson decay, using MCFM g] with the CT10 [5] par- 
ton distribution functions (PDFs). The uncertainty on 
the prediction comes from the PDF uncertainties, eval- 
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uated using the CT10 eigenvector sets, and the QCD 
renormalization and factorization scales, which are var- 
ied simultaneously up or down by a factor of two with 
respect to the nominal value of (mw + ttiz)/2. 

This paper presents measurements of the W^Z pro- 
duction cross-section with the ATLAS detector in pp 
collisions at s/s = 7 TeV. The analysis considers four 
channels of double-leptonic decays W^Z^ i ± vC. + i~ 
involving electrons and muons, i.e. e ± e + e~, /i ± e + e _ , 
e ± /i + /i _ and fj,^ 1 n~ , plus large missing transverse 
momentum. The results are based on an integrated lu- 
minosity of 4.64±0.08fb~ 1 collected in 2011, and super- 
sede the earlier ATLAS results based on a subsample 
of these data [6]. 

The paper is organized as follows: Section [5] briefly 
describes the ATLAS detector and the data sample, in- 
cluding the simulated signal and background samples 
used in this analysis. Section [3] details the definition 
and reconstruction of physically observable objects such 
as particles and jets, and the event selection criteria. 
Section |4] presents the signal acceptance, and Section [5] 
the background estimation. Section[6]presents the mea- 
sured W ± Z production cross-section, constraints on the 
anomalous TGCs, and the fiducial cross-section as a 
function of the Z boson transverse momentum and the 
W^Z diboson invariant mass. 



2 The ATLAS Detector and Data Sample 

The ATLAS detector [7] is a multi-purpose particle 
physics detector operating at one of the beam interac- 
tion points of the LHC0 The innermost part of the de- 
tector is a precision tracking system covering the pseu- 
dorapidity range |?/| < 2.5. It consists of silicon pix- 
els, silicon strips, and straw-tube chambers operating 
in a 2 T axial magnetic field supplied by a supercon- 
ducting solenoid. Outside the solenoid are highly seg- 
mented electromagnetic and hadronic calorimeters cov- 
ering |?7| < 4.9. 

The outermost subsystem is a large muon spectrom- 
eter covering |?/| < 2.7, which reconstructs muon tracks 
and measures their momenta using the azimuthal mag- 
netic field produced by three sets of air-core supercon- 
ducting toroids. This analysis primarily uses the inner 
detector and the electromagnetic calorimeter to recon- 

1 ATLAS uses a right-handed coordinate system with its ori- 
gin at the nominal interaction point (IP) in the centre of 
the detector and the z-axis along the beam pipe. The x-axis 
points from the IP to the centre of the LHC ring, and the 
y-axis points upward. Cylindrical coordinates (r, <j>) are used 
in the transverse plane, <j> being the azimuthal angle around 
the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 9 as rj = — lntan(6*/2). 



struct electrons, the inner detector and the muon spec- 
trometer to reconstruct muons, and the electromagnetic 
and hadronic calorimeters to reconstruct the missing 
momentum transverse to the beam line, E!f lss . The 
E™ lss is corrected to account for muons reconstructed 
by the inner detector and the muon spectrometer. 

W^Z candidate events with multi-lepton final states 
are selected with single-muon or single-electron trig- 
gers. During the 2011 data-taking, the transverse mo- 
mentum (pt) threshold for the single-muon trigger was 
18 GeV. The pt threshold for the single-electron trig- 
ger was initially 20 GeV, and was raised to 22 GeV in 
the latter part of 2011 to cope with increasing instan- 
taneous luminosity. For the W ± Z events that pass all 
selection criteria, the trigger efficiency is in the range of 
(96-99)% depending on the final state being considered. 



2.1 Simulated Event Samples 

Simulated event samples are used to estimate both the 
signal selection efficiency and some of the background 
contributions. The response of the ATLAS detector is 
simulated [3] using Geant4 [S]. 

The production of W ± Z pairs and subsequent de- 
cays are modelled with the MC@NLO [TO1CEI] event 
generator, which incorporates NLO QCD matrix ele- 
ments into the parton shower by interfacing to the Her- 
wig [12] program. The CT10 [5] PDF set is used. The 
underlying event is modelled with the Jimmy [T3] pro- 
gram. 

Background processes for W^Z signal detection are 
jets produced in association with W or Z bosons, 
W + W~ and ZZ pairs, and top-quark production events. 
Alpgen [13] is used to model the W ± /Z + jets and 
Drell-Yan processes for W /Z bosons decaying to e, fi 
and r leptons. Events with multi-jet production from 
heavy-flavour partons are modelled with PythiaB [15] . 
The W + W~ and ZZ processes are modelled with Her- 
wig and Pythia [IB] , respectively. The W ± /Z + 7 
and tt + W ± /Z processes are produced with Mad- 
Graph [17]. The ti and single top-quark events are 
modelled with MC@NLO. Whenever LO event genera- 
tors are used, the cross-sections are corrected to NLO 
or, if available, NNLO matrix element calculations [TUl 

HBH23|. 

Herwig is used to model the hadronization, initial- 
state radiation and QCD final-state radiation (FSR), 
except for the samples generated with Pythia or Mad- 
Graph, for which Pythia is used. Photos [M] is used 
for QED FSR, and Tauola [25J for the r lepton decays. 

Each simulated sample is divided into subsamples 
that reflect the changes in the data-taking conditions 
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in 2011. The average number of interactions per bunch 
crossing, (/i), increased throughout 2011 with the in- 
stantaneous luminosity, and reached a maximum of 17. 
Particles produced in multiple interactions, either co- 
incident with the event of interest or in neighbouring 
bunch crossings, are referred to as 'pile-up' and are in- 
cluded in the simulation. The number of extra interac- 
tions in simulated events is adjusted according to the 
measured (/i) distribution in each data-taking period. 

3 Event Reconstruction and Selection 

The following event selection criteria are applied to the 
events collected with the single-electron or single-muon 
trigger described in Section [2] A primary vertex recon- 
structed from at least three well-reconstructed charged- 
particle tracks, each with pr > 400 MeV, is required in 
order to remove non-collision background and ensure 
good object reconstruction. If an event contains more 
than one primary vertex, the vertex with the largest 
total p\ of the associated tracks is selected. 

3.1 Object Reconstruction and Selection 

Events are selected in the W^Z^ l^vl^l^ channel, 
where the i are either e or /i. The physical objects se- 
lected are electrons, muons, and neutrinos that manifest 
themselves as E T llsa . Contamination from jets, mainly 
due to semileptonic decays of hadrons or due to misiden- 
tification of hadrons as leptons, is suppressed by requir- 
ing the electrons and muons to be isolated from other 
reconstructed objects. 

Muon candidates are identified by matching tracks 
reconstructed in the muon spectrometer to tracks re- 
constructed in the inner detector. The momentum of 
the combined muon track is calculated from the mo- 
menta of the two tracks corrected for the energy loss in 
the calorimeter. To identify muons that traverse fewer 
than two of the three layers of the muon spectrome- 
ter, inner detector tracks that match at least one track 
segment in the muon spectrometer are also included. 
Such muons are referred to as tagged muons to dis- 
tinguish them from the combined muons. The trans- 
verse momentum of the muon must be greater than 
15GeV and the pseudorapidity |?7| < 2.5, using the 
full range of the inner detector. The muon momen- 
tum in simulated events is smeared to account for a 
small difference in resolution between data and simu- 
lation. At the closest approach to the primary vertex, 
the ratio of the transverse impact parameter do to its 
uncertainty (the do significance) must be smaller than 
three, and the longitudinal impact parameter \zq\ must 



be less than 1 mm. These requirements reduce contam- 
ination from heavy flavour decays. Isolated muons are 
selected with a requirement that the scalar sum of the 
Pt of the tracks within AR = 0.3 of the muon, where 
AR = \/(Ari) 2 + (A<j)) 2 , must be less than 15% of the 
muon pt. 



Electron candidates are formed by matching clus- 
ters found in the electromagnetic calorimeter to tracks 
reconstructed in the inner detector . The transverse 
energy Et, calculated from the cluster energy and the 
track direction, must be greater than 15 GeV. The pseu- 
dorapidity of the cluster must be in the ranges |?7| < 
1.37 or 1.52 < \r/\ < 2.47 to ensure good containment 
of the electromagnetic shower in the calorimeters. The 
lateral and transverse shapes of the cluster must be con- 
sistent with those of an electromagnetic shower. The do 
significance must be smaller than 10, and \zo\ must be 
less than 1 mm. To ensure isolation, the total calorime- 
ter Et in a cone of AR = 0.3 around the electron 
candidate, not including the E^ of the candidate it- 
self, must be less than 14% of the electron Et, and the 
scalar sum of the pt of the tracks within AR = 0.3 
of the electron must be less than 13% of the electron 
Pt- The calorimeter response is corrected for the addi- 
tional energy deposited by pile-up. The electron energy 
in simulated events is smeared to account for a small 
difference in resolution between data and simulation. If 
an electron candidate and a muon candidate are found 
within AR = 0.1 of each other, the electron candidate 
is removed. This mainly removes final-state radiation, 
where a photon was misidentified as an electron, and 
also jets from pile-up that were misidentified as elec- 
trons. 



The missing transverse momentum _E™ 1SS is esti- 
mated from reconstructed electrons with \r]\ < 2.47, 
muons with |r/| < 2.7, jets with |7/| < 4.9, as well as 
clusters of energy in the calorimeter not included in re- 
constructed objects with |?7| < 4.5 [27]. The clusters are 
calibrated to the electromagnetic or the hadronic en- 
ergy scale according to cluster topology. The expected 
energy deposit of the muon in the calorimeter is sub- 
tracted. Jets are reconstructed with the anti-fc t jet-find- 
ing algorithm [28 with radius parameter R — 0.4, and 
are calibrated and corrected for detector effects using 
simulation, which has been tuned and validated with 
data. Events that contain jets, with pt > 20 GeV and 
|»7| < 4.9, which are poorly reconstructed as determined 
using calorimeter signal timing and shower shape infor- 
mation, are rejected to improve E™ 133 resolution. 
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3.2 Signal Event Selection 

Events with two leptons of the same flavour and oppo- 
site charge with an invariant mass m« within 10 GeV 
of the Z boson mass are selected. This reduces much of 
the background from multi-jet, top-quark, and W + W~ 
production. Figure HJa) shows the m« distribution of 
the Z candidate in events that pass the complete event 
selection criteria described in this section, except for 
the met requirement. 

Events are then required to have at least three re- 
constructed leptons originating from the same primary 
vertex, two leptons from a Z boson decay and one ad- 
ditional lepton attributed to the decay of a W ± boson. 
To reduce background from Z + jets, the third lepton 
is required to pass more stringent identification criteria 
than required for the leptons attributed to the Z boson. 
The additional criteria imposed on electrons are: a more 
stringent quality requirement for the matched track, a 
requirement on the ratio of the energy measured in the 
calorimeter to the momentum of the matched track, and 
a requirement that transition radiation is detected if the 
candidate traversed the straw-tube chambers. Muons 
attributed to the boson decay are required to be 
reconstructed as combined, and not tagged, muons. Fig- 
ure &b) shows the pr distribution of the third lep- 
tons that pass the additional identification criteria. The 
third lepton is required to have pt > 20 GeV. 

Figure [2jc) shows the E™ lss distribution of the se- 
lected trilepton events. The events have to satisfy E™ lss 
> 25 GeV. 

The transverse mass of the W ± boson is calculated 

as 



(a) 



MX 



2p 



■ T £™ iss (l - cos(/A0)), 



(1) 



where is the transverse momentum of the third lep- 
ton and A(f> is the azimuthal angle between the third 
lepton and the Ef iss . Figure C^d) shows the MXJ dis- 
tribution of the events that reach this stage of the se- 
lection. The observed distribution appears to have 
a narrower peak than predicted by the simulation. The 
events with 70 < Mp < 80 GeV have been scrutinized 
for signs of experimental problems, and no issues were 
found. The limited resolution of the E™ lss measurement 
makes it unlikely that the observed excess is a narrow 
peak. 

Mp is required to be greater than 20 GeV. The 
E^ lss and Mrff requirements suppress most of the re- 
maining background from Z + jets and other diboson 
production. 

In order to ensure that the trigger efficiency is well 
determined, at least one of the muons (electrons) from 
the W ± Z candidate is required to have pr > 20 (25) 
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Fig. 2 (a) Dilepton invariant mass ran °f the Z candidate 
in the events that pass all event selection criteria except for 
the mu cu t. (b) Transverse momentum px of the lepton at- 
tributed to the W candidate, (c) Missing transverse momen- 
tum E™ 133 of the trilepton events, (d) Transverse mass Ms 
of the W candidates. Samples shown in (b), (c), and (d) are 
the candidate events remaining before the cut on the vari- 
able displayed. The stacked histograms are expectations from 
simulation for WZ, ZZ, and Z+7. For Z+jets and ti, the ex- 
pected shape is taken from simulation but the normalization 
is taken from the data-driven estimates. The rightmost bins 
include overflows. 



5 



Table 1 Expected number of signal W^Z^r l^vl+l - events 
after each stage of selection. The first nine rows are computed 
with a simulated W^Z sample scaled to 4.6 fb -1 . Efficiency 
corrections refer to application of correction factors that ac- 
count for the differences in the trigger and reconstruction ef- 
ficiencies between the real and simulated data. The last row 
shows the additional contribution from W^Z — > r ± X. 



eee /iee e/i/i 



Generated 1202 
Muon or electron trigger 1121 
Primary vertex 1118 
Jet cleaning 1116 

Two leptons, ran 219 317 



Three leptons, px 


51.2 


70.6 


74.8 


106.6 


£™ isa > 25 GeV 


40.5 


57.0 


59.2 


86.4 


M™ > 20 GeV 


38.1 


54.1 


55.7 


81.9 


Trigger match 


38.0 


54.0 


55.3 


81.7 


Efficiency corrections 


37.2 


51.8 


54.2 


78.3 


W ± Z — > tX contribution 


1.7 


2.3 


2.4 


3.4 



GeV and to be geometrically matched to a muon (elec- 
tron) reconstructed by the trigger algorithm. These pt 
thresholds arc sufficiently large compared with the trig- 
ger pt thresholds to guarantee that the efficiency is not 
dependent on the px of the lepton. 

4 Signal Acceptance 

The numbers of simulated W^Z events after each stage 
of the event selection, scaled to 4.6 fb _1 , are listed in 
Table[TJ The "Efficiency corrections" row shows the pre- 
dictions corrected for the differences in the trigger and 
reconstruction efficiencies between the measured and 
simulated data. The acceptance increases with the num- 
ber of muons in the final state because the reconstruc- 
tion efficiency for muons is higher than for electrons. 
The additional contribution from W^Z —> t+X, where 
the r decays into an electron or a muon, is shown in the 
last row of Table [TJ 

Table [2] summarizes the systematic uncertainties on 
the expected signal yields. For electrons and muons, 
the reconstruction efficiencies, pt scale and resolution, 
and efficiencies for the isolation and impact-parameter 
requirements are studied using samples of W , Z, and 
J/ip decays. Differences observed between data and sim- 
ulated samples are accounted for, and the uncertainties 
in the correction factors are used to evaluate the sys- 
tematic uncertainties. 

The uncertainties related to lss come mainly from 
the calibration of cluster and jet energies. The effects 
of event pile-up are evaluated from the distribution of 
total transverse energy as a function of (/x) . 

Single-muon and single-electron trigger efficiencies 
are studied in samples of Z — > 11 events. Their effects on 



Table 2 Systematic uncertainties, in %, on the expected sig- 
nal yields. 



Source 


eee 


fiee 




fififi 


fi reconstruction efficiency 


— 


0.3 


0.5 


0.8 


fi pt scale & resolution 


- 


< 0.1 


0.1 


0.1 


fi isolation & impact param. 




0.2 


0.4 


0.6 


e reconstruction emciency 


2.5 


1. ( 


0.8 




e identification emciency 


3.5 


2.3 


1.2 




e isolation & impact param. 


f .5 


1.1 


0.4 




e energy scale 


0.5 


0.3 


0.3 


— 


e energy resolution 


0.1 


0.1 


< 0.1 




J?™ IBS cluster energy scale 


0.4 


0.2 


0.6 


0.2 


E™ 1BB jet energy scale 


0.1 


0.1 


0.1 


0.1 


E™ 13B jet energy resolution 


0.3 


0.3 


0.4 


0.2 


E™ iss pile-up 


0.3 


0.1 


0.3 


0.1 


Muon trigger 




0.1 


0.1 


0.3 


Electron trigger < 


0.1 


< 0.1 


< 0.1 




Event generator 


0.4 


0.4 


0.4 


0.4 


PDF 


1.2 


1.2 


1.2 


1.2 


QCD scale 


0.4 


0.4 


0.4 


0.4 


Luminosity 


1.8 


1.8 


1.8 


1.8 



the W ± Z measurement are small because the presence 
of three leptons provides redundancy for triggering. 

The uncertainty in acceptance due to theoretical 
modelling in the event generator is estimated by com- 
paring MC@NLO with another NLO generator, P0WHEG 
BOX [29] . Uncertainties due to the PDFs are computed 
using the CT10 eigenvectors and the difference between 
the CT10 and MSTW 2008 [3D] PDF sets. Uncertainties 
related to the factorization scale /if and renormaliza- 
tion scale /ii? are estimated by setting fip — fiR and 
varying this value up and down by a factor of two. 



5 Background Estimation 

The major sources of background are summarized in 
Table [31 Data-driven methods are used to estimate the 
background from Z+jets and tt production. Simulation 
is used for the remaining background sources, including 
ZZ, ti+W/Z, and Z + 7 production. Background from 
other sources, such as heavy-flavour multi-jet events, is 
strongly suppressed by the requirement of three leptons 



Table 3 Estimated numbers of background events. The er- 
rors include both statistical and systematic uncertainties. 



Source 


eee 


fiee 


efifi 




Z + jets 


8.8 ±2.8 


3.7 ±2.3 


10.2 ±3.3 


9.1 ±5.5 


ZZ 


3.2 ±0.2 


4.9 ±0.2 


5.0 ±0.1 


7.9 ±0.2 


Z + 7 


1.4 ±0.7 




2.3 ±0.9 




tt 


0.4 ±0.4 


1.7±0.9 


2.3 ±1.1 


2.4 ± 1.2 


ti+ w/z 


0.7±0.1 


1.2 ±0.1 


1.3 ±0.1 


1.6 ±0.1 


Total 


14.5 ±2.9 


11.5 ±2.5 


21.0 ±3.5 


21.0 ±5.6 
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with small do, and is negligible. For studies of anoma- 
lous TGC fSection l6.2l) and of normalized fiducial cross- 
sections f Section 16.31) . the background is estimated sep- 
arately in bins of the transverse momentum of the 
Z boson and the invariant mass mwz of the W ± Z pair. 

For background events with three true leptons from 
vector boson decays, the simulation models the accep- 
tance and efficiency of the selection criteria reliably. 
The main background in this category is ZZ produc- 
tion, in which both Z bosons decay leptonically. The 
background distributions and acceptances are deter- 
mined directly from simulation for this process, and the 
theoretical cross-section is used for normalization. The 
total contribution of the ZZ background is 21.0 ± 0.7 
events, where the error is dominated by the uncertainty 
on the theoretical cross-section. 

Weak boson radiation associated with ti produc- 
tion can also produce three or more leptons and thus 
constitutes a significant background despite its small 
production cross-section. The total contribution of the 
ti+W/Z background is 4.7±0.2 events. Z+j events can 
pass the selection criteria if the photon is misidentified 
as an electron. The contribution from this background 
is estimated from simulation to be 3.7 ± 1.1 events. 

5.1 Z + Jets Background 

Production of a Z boson associated with jets is the 
largest source of background in this measurement. For 
a Z + jets event to pass the event selection criteria, an 
isolated lepton must be reconstructed from one of the 
jets. The extra lepton is usually attributed to the 
boson. 

A lepton-like jet is defined as a jet that passes a few 
basic lepton selection criteria but not necessarily the 
full set of selection (for e) or isolation (for both e and 
n) requirements. An event containing a Z boson and a 
lepton-like jet is a background event if the lepton-like 
jet passes all lepton selection criteria. Those that fail 
the lepton quality or isolation requirements constitute 
a control sample. To ensure that the control sample is as 
similar to the signal as possible, all other event selection 
criteria, including i?™ lss > 25 GeV, are applied. 

In order to estimate the Z + jets background from 
this control sample, the probability / of a lepton-like jet 
passing all lepton selection criteria is estimated in an- 
other control sample: events containing a Z boson and 
a lepton-like jet with £^ iss < 25 GeV. This sample is 
dominated by Z+jets events, and / can thus be directly 
measured. The contributions from other processes are 
subtracted using simulation. Simulation is also used to 
estimate the fraction of the Z+jets background in which 
a lepton-like jet is attributed to the Z boson. 



Table 4 Summary of observed numbers of events iV bs and 
expected signal iV s j g and background TVbkg contributions. 
iVgig includes the contribution from W^Z — > rX. 





eee 


fiee 


efifi 




N obe 


56 


75 


78 


108 


N ais 
AW 


38.9 ±2.1 
14.5 ±2.9 


54.0 ±2.2 
11.5 ±2.5 


56.6 ± 1.7 
21.0 ±3.5 


81.7±2.1 
21.0 ±5.6 



From the combination of the two control samples, 
the total Z+jets background is estimated to be 31.9 ± 
9.2 events. The largest source of uncertainty is the ex- 
trapolation from the E!f lss < 25 GeV sample to the 
^miss > 25 GeV sample, which was studied in simula- 
tion and in dijet data. Also included are the statistical 
uncertainties of the control samples and the uncertain- 
ties on the theoretical cross-sections of the processes 
subtracted from the control samples. 



5.2 tt Background 

A large part of the top-quark background is eliminated 
by the impact-parameter and isolation requirements on 
the leptons, both of which reject lepton candidates orig- 
inating in jets. The rejection factors, however, cannot 
be reliably derived from simulation, and therefore data- 
driven corrections are applied to the simulated ti events 
to estimate them. 

In this analysis, ti events are the only significant 
source of background that does not contain a Z bo- 
son. A control sample enriched in tt background events 
is defined by changing the charge combination of the 
dilepton pair from opposite sign to same sign. All other 
selection criteria are unchanged. Kinematic distribu- 
tions of simulated ti events are similar in shape and 
normalization for same-charge and opposite-charge se- 
lections. The data-to-simulation ratio of the event yield 
in the same-sign sample is 2.2 ± 1.0. This ratio is used 
to scale the ti background predicted in simulation. Us- 
ing this procedure, the total contribution from the ti 
background is estimated to be 6.8 ± 3.2 events. 



6 Results 

The numbers of expected and observed events after ap- 
plying all selection criteria are shown in Table El In 
total, 317 W^Z candidates are observed in data with 
231 ± 8 signal (including final states with r leptons) 
and 68 ± 10 background events expected. There are 206 
W + Z and 111 W~ Z candidates, consistent with the 
expectations of 186 ±11 and 110 ± 6, respectively. 
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Fig. 3 Distributions of the W^Z candidates after all selec- 
tion, (a) Transverse momentum p^, of the Z boson, (b) Invari- 
ant mass rmjvz of the W^Z pair. The shaded bands indicate 
the total statistical and systematic uncertainties of the predic- 
tion. For Z+jets and it, the expected shape is taken from sim- 
ulation but the normalization is taken from the data-driven 
estimates. The rightmost bins include overflow. 



Figure |3] shows the distributions of the transverse 
momentum p^ of the Z boson and the invariant mass 
mwz of the W ± Z pair in selected events. To recon- 
struct mwz, the momenta of the three leptons are com- 
bined with E™ ss , and a boson mass constraint is 
used to solve for the z component p v z of the neutrino 
momentum. This generally leaves a two-fold ambigu- 
ity which is resolved by choosing the solution with the 
smaller \p"\. In 27% of the candidate events, the mea- 
sured transverse mass is larger than the nominal 
mass, and no real solutions exist for p v z . The most likely 
cause is that the measured E^ lss is larger than the ac- 
tual neutrino pt. In this case, the best estimate is ob- 
tained by choosing the real part of the complex solu- 
tions, essentially reducing the magnitude of E™ lss un- 
til a physical solution appears. The data distributions 
are compared with the expected SM signal and back- 
ground. The shapes of the Z + jets and tt distributions 
are taken from simulation and scaled according to the 
data-driven estimates. 



6.1 Cross-Section Measurement 

Two cross-sections are extracted from the number of 
observed events. One is the fiducial W^Z— > l ± vl Jr l~ 
cross-section in a region of final-state phase space de- 
fined by the event selection criteria, the other is the 
total W Z production cross-section. To extract the to- 
tal cross-section, theoretical predictions must be used 
to extrapolate the measured event yield through the ex- 
perimentally inaccessible part of the phase space, intro- 
ducing additional theoretical uncertainties. The fidu- 
cial cross-section is free from such extrapolation, and is 
therefore less sensitive to theoretical uncertainties than 
the total cross-section. 

In order to combine the different channels, a com- 
mon phase space region is defined in which a fiducial 
cross-section is extracted. The common phase space is 
defined as pt^' e > 15 GeV for the leptons from the decay 
of the Z bosons, pt^ e > 20 GeV for the leptons from the 
decay of the W ± bosons, |r?^ e | < 2.5, p v T > 25 GeV, 
\m u - m z \ < 10 GeV, and M]f > 20 GeV, to approxi- 
mate the event selection. In simulated events, the mo- 
menta of photons that are within AR — 0.1 of one of 
the three leptons are added to the lepton momentum. 
In addition, a separation of AR > 0.3 between the two 
leptons of all possible pairings of the three leptons is 
required. This requirement emulates the isolation cri- 
teria applied to the leptons, which tend to reduce the 
signal acceptance for events with very large Z boson 
momenta. The definition of the fiducial phase space is 
identical to that used in Ref. [B] except for the require- 
ment of AR > 0.3 between the leptons. 

For a given channel W^Z— > £^i/£ + £.^, where I is 
cither e or /i, the fiducial cross-section is calculated from 



a wz 



N bs — Vbkg 
J Cdt ■ Cwz 



V r MC 

N MC 



(2) 



where -ZV b s and iVbkg denote the number of observed 
and background events respectively, J Cdt is the inte- 
grated luminosity, and Cwz is the ratio of the num- 
ber of accepted signal events to the number of gen- 
erated events in the fiducial phase space. Corrections 
are applied to Cwz to account for measured differences 
in trigger and reconstruction efficiencies between simu- 
lated and data samples and for the extrapolation to the 
fiducial phase space. The contribution from r lepton 
decays, approximately 4%, is removed from the fidu- 
cial cross-section definition by the term in parentheses, 
where N^ c is the number of accepted simulated W ± Z 
events in which at least one of the bosons decays into 
r, and is the number of accepted simulated W ± Z 
events with decays into any lepton. Since the fiducial 
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phase space is defined by the kinematics of the final- 
state leptons, the calculated cross-section implicitly in- 
cludes the leptonic branching fractions of the and 
Z bosons. 

The total cross-section is defined in the dilepton in- 
variant mass range of 66 < m« < 116 GeV for Z — » 11. 
It is computed as: 



tot 
a WZ 



N 



ba 



bkg 



/ Cdt ■ BwBzAwzCwz 



1 - 



jyMC 



(3) 



where Bw and Bz are the W and Z leptonic branch- 
ing fractions, respectively, and Awz is the ratio of the 
number of events within the fiducial phase space to the 
number of events within 66 < mu < 116 GeV. The ra- 
tio A wz calculated using MC@NLO equals 0.330, 0.332, 
0.333, and 0.338 for the eee, fiee, efXfx, and /i/x/i chan- 
nels, respectively. The differences are due to the FSR 
photons emitted outside the AR = 0.1 cone around the 
electrons. 

Cross-section measurements are extracted using a 
maximum-likelihood method to combine the four chan- 
nels. The likelihood function is defined as 



L(a, x) = n Pois(JV* bs , Nl(a, x) + JVJ(x)) 



(4) 



where Pois(7V* bs , 2V* + N£) is the Poisson probability 
of observing N* hs events in channel i when TV] signal 
and N]; background events are expected. The nuisance 
parameters x affect iV] and N£ as 



k 

j(x)=JVj(0)(l + 5> fc B£), 



(5) 
(6) 



where S k and B k are the relative systematic uncertain- 
ties on the signal and background, respectively, due to 
the fc-th source of systematic uncertainty. 

To find the most probable value of o (fiducial or 
total) the negative log-likelihood function (from Equa- 
tion [4} is minimized simultaneously over a and all the 
nuisance parameters Xk- The final results for the com- 
bined fiducial and total cross-sections are 

0-^ = 92±£(stat.) ± 4(syst.) ± 2(lumi.) fb, 

affi z = 19.0j£f (stat.) ± 0.9(syst.) ± 0.4(lumi.) pb. 

The fiducial cross-section is the sum of the four 
channels. Cross-sections extracted separately for the 
four channels agree within their uncertainties. The un- 
certainties are estimated by taking the difference be- 
tween the cross-section at the minimum of the negative 
log-likelihood function and the cross-section where the 
negative log-likelihood is 0.5 units above the minimum 



Table 5 Systematic uncertainties, in %, on the fiducial and 
total cross-sections. The background uncertainties are split 
into data-driven estimates [Z + jets and tt) and estimates 
from simulation (all other processes). 



Source 


wz 


wz 


\i reconstruction 


0.7 


0.7 


e reconstruction 


2.1 


2.0 


^miss reconstruction 


0.5 


0.5 


Trigger 


0.2 


0.2 


Signal MC statistics 


0.5 


0.5 


Background data-driven 


4.0 


4.0 


Background MC estimates 


0.4 


0.4 


Event generator 




0.4 


PDF 




1.2 


QCD scale 




0.4 


Total 


4.6 


4.8 


Luminosity 


1.8 


1.8 



in the direction of the fit parameter a. The likelihood 
is maximized over the nuisance parameters for each a. 
The systematic uncertainties include all sources except 
luminosity. Correlations between the signal and back- 
ground uncertainties owing to common sources of sys- 
tematics are taken into account in the definition of the 
likelihood. Table [5] summarizes the systematic uncer- 
tainties on the cross-sections from different sources. The 
largest single source of systematic uncertainty is the 
data-driven estimate of the background contributions, 
dominated by that for Z + jets production (±3.8%). 

6.2 Anomalous Triple Gauge Couplings 

General expressions for the effective Lagrangian for the 
WWZ vertex can be found in Refs. |3 1U32) . Retaining 
only terms that separately conserve charge conjugation 
C and parity P, the Lagrangian reduces to 

Cwwz = i g?(Wl u W»Z v - W^W^Z") 
gwwz [ 



+ K Z WlW v Z» v + 



A; 



rr pfj, r y V 



where gwwz = — ecot 9w, e is the elementary charge, 
Bw is the weak mixing angle, and Z^ are the W 
and Z boson wave functions, X^ v = d^Xy — dyX^ for 
X = W or Z, and gf , kz, and Xz are dimensionlcss 
coupling constants. The SM predicts gf — 1, Kz — 
1, and Xz = 0. This analysis sets limits on possible 
deviations of these parameters from their SM values, 
i.e. on Agi = gf — 1, Akz = «z — 1, and Xz, known as 
the anomalous TGC parameters. The W ± Z production 
cross-section is a bilinear function of these anomalous 
TGCs. 

To avoid tree-level unitarity violation, the anoma- 
lous couplings must vanish as s, the four-momentum 
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Table 6 Expected and observed 95% confidence intervals on 
the anomalous couplings Ag^ , Akz, and Xz- The expected 
intervals assume the Standard Model values for the couplings. 



2000 
P z [GeV] 

Fig. 4 Transverse momentum p^, of the Z boson in W^Z can- 
didate events. Data are shown together with expected back- 
ground and signal events, assuming the Standard Model. Ex- 
pected events in the case of anomalous TGC without form 
factor are also shown for illustration. The last bin is short- 
ened for display purposes. 



squared of the W ± Z system, approaches infinity. To 
achieve this, an arbitrary form factor may be intro- 
duced [35]. Here the dipole form factor adopted is 



a{s) 



«o 



(1 + s/A 2 ) 2 



(8) 



where a stands for Agf , Akz, or Xz, ao is the value 
of the anomalous coupling at low energy, and A is the 
cut-off scale, the scale at which new physics enters. The 
results are reported both with and without this form 
factor. 

Since an enhancement in the cross-section due to an 
anomalous coupling would grow with s, measurement 
sensitivity to anomalous TGCs is enhanced by binning 
the data in a kinematic variable related to s. The trans- 
verse momentum of the Z boson provides a natu- 
ral choice for such binning as it is strongly correlated 
with s and can be directly reconstructed from the mea- 
sured lepton momenta with good precision. The data 
are therefore divided into six bins in p^ of width 30 GeV 
followed by a wide bin that includes 180-2000 GeV. 

MC@NLO [IT] is used to generate W ± Z events with 
non-SM TGC. The generator computes, for each event, 
a set of weights that can be used to reweight the full 
sample to any chosen set of anomalous couplings. This 
functionality is used to express the predicted signal 
yields in each bin of p^ as a function of the anoma- 
lous couplings. Figure S] shows the p^ distribution of 
the selected events together with the SM prediction. 
Also shown for illustration are predictions with non- 
zero anomalous couplings without form factor: each cou- 
pling is increased to the expected 99% confidence-level 
upper limit while keeping the other two couplings at 
the SM value. For this plot the 99%, rather than 95%, 





Observed 


Observed 


Expected 




A = 2 TeV 


no form factor 


no form factor 




[-0.074,0.133] 


[-0.057, 0.093] 


[-0.046,0.080] 


Ak z 


[-0.42,0.69] 


[-0.37,0.57] 


[-0.33,0.47] 


Az 


[-0.064,0.066] 


[-0.046, 0.047] 


[-0.041,0.040] 



confidence-level upper limits are used to accentuate dif- 
ferences in shape. As expected, the largest deviations 
from the SM are in the last bin of p^,, while the devia- 
tions in the lower-p^ bins depend on which coupling is 
varied. 

Frequentist confidence intervals are obtained on the 
anomalous couplings by forming a profile likelihood test 
incorporating the observed number of candidate events 
in each p^ bin, the expected signal as a function of 
the anomalous couplings and the estimated number of 
background events [33] . The systematic uncertainties 
are included in the likelihood function as nuisance pa- 
rameters with correlated Gaussian constraints. A point 
in the anomalous TGC space is accepted (rejected) at 
the 95% confidence level if less (more) than 95% of ran- 
domly generated pseudo-experiments exhibit a value of 
the profile likelihood ratio larger than that observed in 
data. 

Table [5] summarizes the observed 95% confidence 
intervals on the anomalous couplings Agf , Akz, and 
Xz, with the cut-off scale A = 2 TeV and without the 
form factor. The limits on each anomalous TGC param- 
eter are obtained with the other two anomalous TGC 
parameters set to zero. The expected intervals in Ta- 
bleware medians of the 95% confidence-level upper and 
lower limits obtained in pseudo-experiments that as- 
sume the SM coupling. The widths of the expected and 
observed confidence intervals are dominated by statisti- 
cal uncertainty. Figure [S] compares the observed limits 
with the Tevatron results [Mll3"5"] . 

The 95% confidence regions are shown as contours 
on the {Agf , Akz), {Agf,Xz), and (Akz,Xz) planes 
in Figure El In each plot the remaining parameter is set 
to the SM value. The limits were derived with no form 
factor. 



6.3 Normalized Fiducial Cross-Sections 

The effective Lagrangian adopted in the TGC analysis 
in Section 16.21 allows us to probe non-SM physics with 
little model dependence. An alternative approach is to 
measure kinematic distributions, such as the p^* spec- 
trum, that could be compared with model-dependent 
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Fig. 5 95% confidence intervals for anomalous TGCs from 
ATLAS (this work), CDF [31], and DO [35]. Integrated lumi- 
nosity, centre-of-mass energy and cut-off A for each experi- 
ment are shown. 



theoretical predictions. For this purpose, it is neces- 
sary to convert the measured distributions to the un- 
derlying true distributions by unfolding the effects of 
the experimental acceptance and resolution. The iter- 
ative Bayesian unfolding proposed by D'Agostini [55] 
is applied here. An implementation of this technique 
has previously been used by ATLAS to unfold the pr 
spectrum of inclusively produced W bosons [37] ■ 

In the unfolding of binned data, effects of the ex- 
perimental acceptance and resolution are expressed in a 
response matrix, each element of which is the probabil- 
ity of an event in the i-th true bin being reconstructed 
in the j-th measured bin. In iterative Bayesian unfold- 
ing, the response matrix is combined with the measured 
spectrum to form a likelihood, which is then multiplied 
by a prior distribution to produce the posterior proba- 
bility of the true spectrum. The SM prediction is used 



as the initial prior, and once the posterior probability 
is obtained, it is used as the prior for the next iteration 
after smoothing. The spectrum becomes insensitive to 
the initial prior after a few iterations. The number of 
iterations is adjusted to control the degree of regulariza- 
tion [36] . The differences between successive iterations 
can be used to estimate the stability of the unfolding 
method. 

To achieve stable unfolding, that is, without exces- 
sive sensitivity to statistical fluctuations in data or to 
details of the unfolding technique, the measured quan- 
tity must be a good approximation to the underlying 
true quantity: the response matrix must be close to di- 
agonal. The distribution used in the TGC analysis 
is a natural choice that has good resolution and sensi- 
tivity to possible new physics. The fractions of W^Z 
events that migrate between two p^ bins are 2-7%. 

In addition to p^, the distribution of the diboson 
invariant mass mwz is also measured. The resolution 
of the reconstructed mwz is limited by the E™ lss res- 
olution. To avoid large bin-to-bin migration and achieve 
stable unfolding, three mwz bins are used: 170-270 GeV, 
270-405 GeV, and 405-2500 GeV. With this binning, 
the fractions of events that migrate between two mwz 
bins are 13-17%. 

Figure [7J shows the fiducial cross-sections extracted 
in bins of p^ and mwz, normalized by the sum of all 
bins. Comparison with the SM prediction shows good 
agreement. The corresponding numerical values are pre- 
sented in Tables [7] and [5] for and mwz, respectively. 

The dominant source of uncertainties on the normal- 
ized cross-sections is statistical. The statistical uncer- 
tainties are determined by a Monte Carlo method. Two 
thousand pseudo-experimental spectra are generated by 
fluctuating the content of each bin according to a Pois- 
son distribution. The unfolding procedure is applied to 
each pseudo-experiment, and the r.m.s. of the results is 
taken as the statistical uncertainty. The systematic un- 
certainties are evaluated by varying the response matrix 



11 



Table 7 Normalized fiducial cross-sections and uncertainties in bins of p 2 , . 



P? [GeV] 


[0, 30] 


[30, 60] 


[60, 90] 


[90, 120] 


[120, 150] 


[150, 180] 


[180, 2000] 


Uncertainty 


0.231 
0.034 


0.350 
0.039 


0.230 
0.033 


0.065 
0.019 


0.045 
0.015 


0.042 
0.014 


0.038 
0.013 



(a) 
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Fig. 7 Normalized fiducial cross-sections Acr^ z /a f ^ z in bins 
of (a) p 2 , and (b) mwz compared with the SM prediction. 
The total uncertainty contains statistical and systematic un- 
certainties added in quadrature. 



for each source of the uncertainty, and combining the 
resulting changes in the unfolded spectrum. Because of 
the normalization, the results are affected only by the 
uncertainties that depend on or mwz ■ The stability 
of the unfolding procedure is tested in two ways: firstly 
by comparing the unfolded spectra after two and af- 
ter three iterations, and secondly by checking that the 
true variable distribution is correctly reproduced from 
a simulated sample generated with non-zero anomalous 
couplings. 



Table 8 Normalized fiducial cross-sections and uncertainties 
in bins of mwz- 



m wz [GeV] 


[170, 270] 


[270, 405] 


[405, 2500] 


Uncertainty 


0.568 
0.038 


0.283 
0.030 


0.149 
0.027 



data sample with an integrated luminosity of 4.6 fb , 
collected with the ATLAS detector at the LHC. The 
candidate W^Z events were selected in the fully lep- 
tonic final states with electrons, muons, and large miss- 
ing transverse momentum. In total, 317 candidates were 
observed with a background expectation of 68 ± 10 
events. The fiducial and total cross-sections are deter- 
mined to be 

= 92±£(stat.) ± 4(syst.) ± 2(lumi.) fb, 



and 



tot 

a wz 



19.01]; ICstat.) ± 0.9(syst.) ± 0.4(lumi.) pb, 



respectively where the fiducial cross-section is defined 
by Px' e > 15 GeV for the leptons from the decay of 
the Z bosons, p^ e > 20 GeV for the leptons from the 
decay of the W ± bosons, |^' e | < 2.5, p£ > 25 GeV, 
\m u - m z | < 10 GeV, Mjf > 20 GeV, and AR > 0.3 
between the two leptons of all possible pairings of the 
three leptons. These results are significantly more pre- 
cise than the earlier ATLAS measurement [5] which 
this paper supersedes. The total cross-section is con- 
sistent with the SM expectation of 17.613'J pb. Limits 
on anomalous triple gauge couplings have been derived 
based on the observed distribution. The 95% confi- 
dence intervals are 

Agf <= [-0.057,0.093] 
Ak z E [-0.37,0.57] 
X z E [-0.046,0.047] 



7 Conclusion 

Measurements of W^Z production in proton-proton 
collisions at y/s = 7 TeV have been presented using a 



without a form factor. The limits are again more strin- 
gent than the earlier ATLAS measurement. Normalized 
fiducial cross-sections have also been presented in bins 
of and mwz, and are in good agreement with SM 
predictions. 
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